TAUP-2892/09 

Holographic Dual of QCD from 
Black D5 Branes. 



^Benjamin A. Burrington, * Jacob Sonnenschein 

^Department of Physics, 
University of Toronto, 
Toronto, Ontario, Canada M5S 1A7. 

* School of Physics and Astronomy, 
The Raymond and Beverly Sackler Faculty of Exact Sciences, 
Tel Aviv University, Ramat Aviv, 69978, Israel. 

Abstract 

We study the dynamics of probe D7 D7 flavor probe branes in the back- 
ground of near extremal D5 branes. This model is a holographic dual to a 
gauge theory with spontaneous breaking of a U{Nj)l x U {Nf)ii chiral symme- 
try. The spectrum of two such D7 D7 embeddings, contains a single massive 
4D meson coming from the world volume U{\) gauge field, the pion, and a 
single massive 4D scalar meson coming from fluctuations of the embedding of 
the brane. In addition, there are continuum five dimensional states due to the 
finite height of the effective potential in the radial direction. We investigate 
baryons in this model, and find that the size is stabilized due to the Chern- 
Simons term in the D7 world volume action. The model admits a Hagedorn 
temperature of 2^ where R is the radius parameter in the D5 branes metric. 
We investigate the pattern of chiral symmetry breaking in the deconfined phase 
as a function of the asymptotic separation of the branes L. We find that for 
< L ^ 1.068i? that chiral symmetry is restored, and that chiral symmetry 
is broken for L outside this window. We further argue that the solutions with 
L < ^ are only classically stable, and in fact no D7 embedding exists with 
these boundary conditions. 



1 Introduction 



QCD at strong coupling has several properties of interest: Confinement, the exis- 
tence and spectrum of bound states, and chiral symmetry breaking to name a few. 
Studying these phenomena from the standpoint of quantum field theory is difficult, 
and impossible when confined to the realm of perturbation theory. 

Recently, much progress has been made by constructing brane background models 
with a tunable parameter, such that in one limit they become pure (weakly coupled) 
SU{Nc) QCD, and in the opposite they become a weakly coupled string theory [1]. 
Such constructions allow one to study a model that is in the same universality class as 
pure YM theory via gauge gravity duality [2] (for a review see [3] ) . Karch and Katz [4] 
introduced a holographic dual of dynamical flavored fundamental quarks by adding 
flavor probe D-branes into the gravity background. By taking the number of flavor 
branes to be much smaller than Nc, Nf <^ Nc, the backreaction on the background 
is avoided. This is the analog of the quenched approximation in the gauge theory. 
This was done originally for non-confining backgrounds. In [5] a similar scenario 
was proposed in a confining setupll] The incorporation of chiral symmetry and its 
spontaneous breakdown has been done in the seminal papers of Sakai and Sugimoto 
(SS) [7,8]. Their model which has been extensively studied in the literature. 

In the SS model, the background is taken to be a near extremal D4 with one of 
the spatial coordinates compactified, which we will call X4, and its radius of com- 
pactification is Rx^- For energies E -C one expects 4D dynamics contaminated 

-n-x4 

with KK modes. In the low temperature regime, the correct place to put the black- 
ening factor is in front of the dxl component of the metric, giving the u,X4 plane 
the geometry of a (semi-infinite) "cigar". The probe flavor branes are taken to be 
(Nf) D8 branes, and are parallel to some combination of the radial direction, pa- 
rameterized by u, and X4. In this model, the realization of chiral symmetry breaking 
U{Nf)fi X U{Nf)L — * U{Nf) is geometric, and easy to visualize: the stack of D8 
branes must connect smoothly at the bottom of the cigar, and so connect smoothly 
out to D8 branes at infinity, making a U-shaped profile in the u, X4 plane. At large 
u, the UV, one sees a disconnected D8 D8 pair with U{Nf)fj x U{Nf)L which is 
broken at small u, the IR, where one sees only the diagonal U{Nf) symmetry re- 
maining. A variety of physical properties of meson and baryon physics has been 
extracted from the model. These include the massive meson spectrum, the massless 
Goldstone pions [7], certain decay rates [8], thermal behavior of hadrons [9], [10], as 
well as describing certain condensation phenomena in the theory [11-17]. It is clear, 
however, that the SS model is far from being the ultimate gravity dual of QCD. It 

^ For references of additional papers that discussed holographic dynamical quarks see [6] 
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suffers from several drawbacks like the fact that its UV completion is not known, 
like discrepancies between its hadronic spectrum and that observed in nature. Thus, 
one is instructed to further construct holographic laboratories. The present paper 
follows these lines. 

Here, we will instead take a near extremal D5 background with two of the world 
volume coordinates compactified (x4,X5), and take the flavor branes to be probe D7 
branes. We take probe D7 branes will be parallel to some combination of u,X4^,X5] 
some aspects of this setup have been investigated in [18], and a close analog inves- 
tigated in [19]. We expect that the low energy dynamics to be once again 4D, and 
given the arguments of [1] (and applied to the D5 background in [19]), should also 
be in the same universality class as pure YM theory, excluding the flavor branes. 
Including the flavor branes adds chiral fermions to the low energy weak coupling 
description, and so the holographic setup is dual to a theory in the same universality 
class as pure large Nc QCD. 

There are, however, some key differences between extremal D4 and D5 branes. 
Most distinctly is the fact that the decoupling limit does not raise a large barrier, 
decoupling the bulk modes. Rather, the radial effective potential becomes infinitely 
long, but of finite hight. We find that this is true for the effective potential of the D7 
brane fluctuations as well. This gives that the effective potential has finitely many 
bound states, and so unfortunately only a finite number of 4D mesons. We will, 
however, see the same type of geometric understanding of chiral symmetry breaking 
as described above. 

To begin, we will now briefly review the near horizon limit of the near extremal 
D5 background we wish to study. The metric is given by 

u R dii^ 

= t; iVuudx^'dx'' + dxl + dxifiu] u^)) H + Rud^l (1) 

R u f{u;uA) 

where dQ^ is the metric of the unit 3 sphere, and 

/(«;«a) = (i-^). (2) 

Furthermore, there is a non trivial dilaton and 3 form flux 

exp(j) = gs — , Fs = iZg (3) 

R 9s 

where VL^ is the volume form of the unit 3 sphere. The parameter R is related to 
string parameters in the usual way R^ = QsNcQ'. We will be considering the case 
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where both and are compact 



X4 = X4 + 27rR^^, X5 = X5 + 27ri?a.g. (4) 

The parameter R^^^ is in fact not independent of the other scales aheady mentioned, 
given that we want the solution to be smooth around u = Up^. One can see this most 
easily by changing to variables 

= u\ + z^. (5) 

In these variables, the metric reads 



ds^ = ^ (7]^,dx''dx'' + dxl) + RJul + z^dnj + \ ^ ^ (6) 



and clearly, for the metric to be non singular around z = 0, we need to take 

= X5 + 2'kR, i.e. Rx^ = R. (7) 

Thus we see that unlike the SS here ua is a free parameter not related to R. In fact 
this is the situation classically. Based on similar situations [20], it is plausible that 
quantum mechanically the system will not be stable for any ua. However, for the 
case discussed in section 3 this instability does ont occur. One should also note that 
the above background is confining, as ^JguQxx is finite. 

For this background, we read off the relation between field theory and string 
theory quantities 

2 /o n3 / 2 '^T^QsOi' 



% = (27r) Qsa', 



RR 



x4 



1 Ua 

Tst = ^^/9u9^\u=u^ = (8) 

The fact that is a free parameter implies that so is the string tension. Further, we 
should comment that this background is S-dual to the near extremal NS5 brane. Such 
NS5 brane constructions have been argued to be the holographic dual of little string 
theories (LSTs) [21] (for a review, see [22]), however in the case at hand we are in the 
opposite range of validity {qs <^ 1). Further, because S duality works non trivially 
on the D7 brane probes, our theory will not have a straightforward UV completion. 
The use of holographic techniques for the NS backgrounds, [23-25], should, however, 
be easily be adapted to our scenario (and in some cases be identical). 
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We organize the remainder of the paper as follows. In section 2 we discuss the 
general behavior for D7 embeddings in the confined phase, and briefly discuss the 
eigenvalue problem to describe the masses of mesons in this general case. In section 
3, we consider D7 embeddings into the extremal limit of the background, which 
is also sufficient to describe D7 embeddings that do not come near to the end of 
the radial coordinate at u = u\. We find that these embeddings obey a simple 
scaling behavior, which we describe generally enough to be applicable outside of 
these particular models. We also consider fluctuations of the embedding functions 
and the world volume gauge field to describe the spectrum of mesons in this limit. 
In section 4, we describe an embedding that falls all the way to m = ma, and again 
consider fluctuations of the embedding functions and world volume gauge field. We 
also consider baryons in this scenario. Finally, in section 5 we consider the deconfined 
case, and consider the pattern of chiral symmetry breaking as a function of the 
boundary conditions. In section 6 we conclude, and suggest some future lines of 
investigation. 



2 D7 embeddings 

In this section we will consider embeddings of D7 branes into the above background. 
Qualitatively, we will find embedding of the form shown in figure [H This background 




Figure 1: The configuration of D7 branes in the near extremal D5 brane background, 
following the above considerations. 

is similar to the paperclip or hairpin background (see for instance [26] and [27]). The 
embeddings we consider are transverse to some combination of the directions, 
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and so the Chern Simons (CS) term will be unimportant in determining the solutions 
of the equations of motion, and we are left with the Dirac-Born-Infeld (DBI) action 




Sd7 = -K7 / d«eexp(-0)^/ ( -det l^g,,^^^jj (9) 

where we have used the notation to represent world volume coordinates of the D7 
brane, and is the D7 tension. We use an Ansatz of the form 

x\^) = « = yu, i = sphere coordinate (10) 
X5iO = x^ip), XiiO = X5(p), u(0 = u{p). 

In the above notation, we have picked a "static gauge" for the and sphere direc- 
tions. The final world volume coordinate, p, parameterizes a path in X4^,X5,u space. 
One can show directly that the full non linear equations of motion for the and 
sphere directions are satisfied for the above Ansatz. Using the reparameterization 
invariance in p, we make one final choice 

uip) = p. (11) 



We now reduce the action to 1 dimension (p) and find the relevant equations of 

d_ 

dp 



motion. In the following we will abbreviate = d. Using this notation, the action 



becomes 



Sm = / dpp\ {dx,Y + (9x5)V(p;^a) + r- (12) 

The equations of motion are quite simple, and follow from the Noether charges 
associated with translation invariance in X4 and X5. These read 

= P,< (13) 



(8x4) + (8x5) /(p;ua) 



p^/(p;"a) 



P'f{p;u,)8x, ^^^^3_ (14) 



^ {8X4)^ + {8X5)^ f{p] Ua) + ^ 



/(p;"a) 
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This is a completely integrable system: one may simply solve for 8x4 and 8x5 and 
then integrate with respect to p. Doing so, we find 



X4 



dpP^ 



UA 



^] -PI 



Pi 



/(p;"a) 



(15) 



X5 



dpP^ 



\ 



p2 Pj 
4 /{p;nA) 



(16) 



The above expressions represent only one branch of the solution: in the end one must 
match an identical branch connecting at the minimum value of uq = po- One may 
check explicitly that the field equations of motion arising from action (Q are met for 
the solution given by (ITS!) . (fT6|) and (ITT!) . 

There is a similarity to the solution for x^ and the solution given in [9], except 
that a different power of p appears with the (also different) function /. In fact, for 
the case P4 = 0, we find that the minimum value of u is in fact given by uq, with 

— ^ \//(mo; wa) which is a very similar expression to that given in [9] . 

The new piece of information is given by P4, and this changes the situation in an 
interesting way. First, we note that the "turn around point" for the embedding is 
given where dx^ or dx^ blow up (i.e. m = p is not varying while both X4 and X5 vary 
allot). Note that for both x^ and X5 this happens at the same place: either where 

f{p]UA) ^ or where ( f — ) —Pi — jr^] ~^ 0- The question then becomes. 



\ \MA / 4 f(p;uA) I 

which one goes to zero first? 

To help address these questions, we first change to variables 

P 

P = —■ 

Ua 

Under this transformation we are considering where 



(17) 



/(p;nA) = /(p;l) = 



P 



p2 _ 
-'4 



p2 



fip; 1) 



0. 



(18) 



We can now see which one goes to first: it is clearly the second line (when P5 7^ 0). 
Imagine coming in from very large positive p and decreasing this value. The first 
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line goes to zero only when p = 1. However, in such a limit, the second line already 
has a zero. This is because of the / in the denominator of the last term is going to 
zero (from the positive side), and so at some point this will compensate for the value 
of — Pi, and have a zero at a value of p > 1. 

The only special case is when P5 = 0. When this is the case, one of the functions 

goes to zero at max ^P/ , 1 j . Note that this has the peculiar feature that the end- 
points of the D7 brane (a function of P4) may not be at the same point, but none 
the less, the brane falls all the way to the bottom of the cigar u = ma- This ceases 
once IP4I > 1. If IP4I > 1, then the embedding isn't antipodal any more in X5. In the 
Xs, u plane, the embedding never reaches the bottom point, and so then never needs 
to be matched to a solution going through the origin (we will see this more clearly 
we plot the asymptotic separation distance as a function of P4 and P5). 

One can check the above discussion of uq = po = poU\ simply by plotting the 
results. We define po by the equation 

Clearly, po is a function po(P4, P5) and so we may plot this (we plot the largest real 




Figure 2: A plot of the largest real root, pq^P^^.P^) = no(P4,P5). Note all values 
are larger than 1. The wire frame in the right hand plot corresponds to the large Pj 
asymptote po(P4, P5) ^ (P4 + ^1)^- 

root), and compare it with the value 1 in figure O Further, we can see the special 
case of P5 = arising, and the flattening of the function that happens along the line 
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interval IP4I < 1: indeed max yP^ , Ij is the limiting value of the solution along the 

line P5 = 0. Further, for large values of P4 and P5, the value of the zero po grows as 
well, and so the / in the last part of the expression can be approximated by 1. This 
gives that the approximate value of the zero is given by po = (P4 + Pi)'^' ■ We plot 

this limiting value and show that it does indeed asymptote (see figure [2]). Finally, 

1 

one can show that this surface exactly corresponds to the curve max(P/ , 1) along 
the line P5 = by plotting a side view (see figure [3]). 




Figure 3: The left hand graph is a side view of the plot of the root of the above 

rational function, plotted as a function of P4, P5. The black wire frame plot is the 

J 1 

asymptote po = {Pi + -^5^)^ = Pi when P5 = 0. The right hand graph is a plot of 
Po(P4, P5) zoomed in on the region P5 = 0, IP4I < 1. One may see that this surface is 
smooth on the line interval P5 = 0, IP4I < 1 except the endpoints P5 = 0, P4 = ±1. 

One may be curious how po ^ 1 as P5 — on the line interval IP4I < 1. Clearly, 
it must relax as P|, and so must be smooth around P5 = 0. We plot a zoomed in 
version of the plot of po and find that this is indeed the case. The only singular 
points on the surface (with conical singularities) are at (P5, P4) = (0, ±1). See figure 

m 

Now that we have found (qualitatively) the minimum value to which the brane 
falls, we also wish to relate the parameters P4 and P5 to brane separation at u = 00 
(in the coordinates and 0:5). We will call the asymptotic separation in Ax 4 = L4 
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and Axs = L^. In the end it will be these separations that we wish to consider the 
"boundary conditions" imposed on the branes, so that the fluctuations around these 
embeddings will satisfy Dirichlet boundary conditions. Given the above expressions, 
we see that 



where the factor of 2 comes from "gluing" the second branch of the solution on. 
Again, finding analytic solutions to the above is problematic, but we can graph the 
results function of P4 and P5, see figure |H 

We first consider the graph of L^/R. Before making this graph, one may think 
that we need to be careful about finding the actual separation distance L5: one 
could have a solution that connects through smoothly the point u = u^, and so the 
second branch actually lies displaced by ttRx^ = ttR. However, the integral for L^^/R 
appears to have this built in. The only cases where the turn around point of the 
branes lies at Uq = uj>^ is on the interval — 1 < P4 < 1 and P5 = 0, and so we focus 
on this region. While we do not offer an analytic proof, we may zoom in on this 
part of the graph, and be convinced that L5/P — > ivr (see figure [5]). One may also 
be convinced simply by noting that if P4 = P5 = 0, then the second branch must 
be glued on in antipodally after reaching u = ua to avoid a discontinuity in the 
tangent vector in m,X5,X5 space. If this is the case, we see that the integrand in the 
expression for L5 must collapse to a delta function with coefficient ivr in the limit 
that P5 ±0. For generic values of p the limit P5 — >■ sets the integrand to zero, 
however at p = po it blows up to infinity with the integral underneath remaining 
constant. This is sufficient to show that the integrand goes to a delta function in p 
for the limit P5 ^ 0. This also explains why one misses this if one passes the limit 
P5 = under the integration: the limit implied by the integral and the limit P5 ^ 
do not commute, and so one must put in the additional vr of separation by hand. 
Further, because the values ivr are identified, the surface is actually continuous. 

Next we consider the graph for L4. In figure H] there is a divergence in L4 around 
P5 = and P4 ±1. This is easy to understand. If one looks at P5 = and 
P4 = 1, the integrand goes as around p = po = 1, and so the integral is log 
divergent. Because L4 blows up around P5 = 0, P4 = ±1 we may wish to exclude 
multiple wrappings of the D7 around X4. We further restrict the total X4 winding to 




(20) 




(21) 
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5.0 5.0 .5.0 " -5.0 




Figure 4: Plots of L4^{Pi, P5)/R and L^i^P^, P5)/R (labeled simply L4 and L5 in the 
graphs). The view on the left gives the same perspective as earlier graphs, while 
the view on the right is slightly easier to see what is happening near the surface 
L4(P4, P5)/i? = and L^iP^, P^)/R = 0. 

be n/Rx^ (half the circumference), knowing that there is another solution where the 
same boundary conditions are satisfied, and less distance has been traversed by the 
brane. Hence, one may wish to exclude a small region around P5 = and P4 = ±1 
whose boundary is given by IL4I/P = nRx^/R. 

The above considerations are what lead to the diagram in figure [H 
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-2 



-1.0 



P5 



Figure 5: Plots of L^i^P^, P^)/ R zoomed in around the interval. Note that in the side 
view the maximum and minimum appears to approach ivr. 

2.1 Mesons in the model. 

The mesons of this model are described by strings ending on the D7 branes [28]. 
The high spin mesons are given by semiclassical strings which end on the D7 at the 
minimal value of the brane embedding mq, fall towards the minimum value of ua and 
lay flat at this value and finally climb back up to meet the D7 again at uq. The low 
spin (vector) mesons are given by fluctuations of the D7 U{N) gauge field strength, 
and the scalar mesons are given by the fluctuations of the embedding functions. 
To compute any of these fluctuations, we will need the puUback metric 



dSp 



P_ 
R 



■q^j^ydx^dx" + ( dxl + dxlf{p; ua) + 



R^ 1 

/(p;ma) 




(22) 



where we have defined the useful function 



P 



,6 




(23) 



For later use, we will also need 




(24) 
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and we find it convenient to use the notation 



With such an assignment, the solutions above read as 

1 HP 
dx, = RP,^= _ = ^^) (26) 

/pV(p; 1) (p^ - PI - 7§y^ ^ 



1 /? p 

fe5 - flA , ^ = JTji^,VW)- (27) 



4 /(P;l) 



2.1.1 Vector mesons 



First, we compute the effective action for the fluctuations associated with the gauge 
field Fab which we take to be only a function of /), and further that components 
of F with indices along the sphere directions are zero. With this, we take the DBI 
and expand it to second order in F |^ 



Sd7 = -^7 j d^^exp(-0)y^- det {g{p)ab + Fab) 
= -1^7 j d^^exp(-0)y^- det {g{p)ab) 



1 + g^^'^Fa^b, + \Fa,b.Fa,b, [gi''' g''^'^ ' gi'"' g?'') + ■ ■ ■ ) (28) 

and now because F is antisymmetric, while g is symmetric, we find that to leading 
order 



-K-j j J dx^dp^p^y^ (^1 + ^Fa,b,Fa,b,g;"''g:^ 



f)2 

P 



.(29) 



■^We need not worry about the CS term because the puUback of Ae (generated by the D5 
background) onto the D7 worldvolume is zero. Note also that we are takmg F to be a U(l) gauge 
field. One may generalize this discussion to a U(N) field by promoting _F to a matrix, and taking 
a trace. 
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Hence, the quadratic action we wish to consider is 



Sf2 = — ^J^^ f dx^dpu\Rp^ \/ '^{p) 
^9s J 



E 2 



p Ma 



Ma7(p) 



4^. 



and in the above we have used rj^^ to raise and lower indices. 
The equations of motion following from the above action are 



d^F^'p = 



d^F^f, = 0. 



We expand the gauge field in the following way 



(30) 



(31) 
(32) 

(33) 



We further gauge fix d^B^jf^ = 0. Under these conditions, and for the above gauge 



(34) 



choice, we see that the equation fl32|) becomes 



where m„ is defined by d^d^By = m^By. We see that 0„ = by this equation of 
motion except for the case where m„ = 0. We choose to index such a mode mo = 0. 
vTo is the massless pion. 

We now note that the remaining equation of motion (13T]) reads 



E 



d 



0. 



(35) 



One can show that there are no massless normalizable modes of the vectors By, and 
so we find that 



d 



^ 2 

^'0 oc — . 
p2 



Thus, we find that the eigenvalue problem for the ipn reduces to 



-d ( P^d^.^ 

^ 2 



(36) 



(37) 
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where the 4D mass is given by •mA,n = ■ 

Unfortunately, the eigenvalue problem for the fluctuations of the embedding co- 
ordinates are much more involved because we do not have a good set of orthogonal 
coordinates to the brane. However, we will study these in two special cases below. 



3 Large uq ^ u\ limit. 

In this section, we consider the limit where uq ^ ua. This limit is the same as 
the extremal limit of the background because f{u; n^) can be approximated by 
f{u\utC) = 1- As mentioned in the introduction this case makes sense also quan- 
tum mechanically. In fact it is similar to the cases studied in [29]. In such a limit 
the periodic identification of and now both become arbitrary, and further there 
is an additional 5*0 (2) symmetry rotating the X4,,x^ plane, which the solutions will 
transform under in a simple way. To consider the limit mq ^ ma, we cannot take 
a strict ma ^ limit of equations (fT3l) and (fT^ because this would scale the right 
hand side of these equations to zero (when they should be constants). Therefore, we 
absorb the u\ into Pj, via pi = u\Pi (the parameters pi have units). One may think 
of this as a scaling where ma — ^ while Pj — > oo with = u\Pi remaining fixed. 
This is the connected to the fact that for large Pj the solution is in a regime where 
no > ua. 

We consider equations (fT3ll and (fill) with the above substitutions and set ma = 
and find 



p^dx4 



^{dx,f + {dx,f + f 
p^dx^ 



{dx,Y + {dx,Y + f 



As before, we solve for dx^ and 9x5, and find 



= P4 (38) 
p,. (39) 



^ Rpi arcsin 

X4 = / dpRpi^l ^ — ^ = , ^ ^ + X4,o (40) 

V ^4 + 




,3 



^ Rp5 arcsin i - 

X5= / dpRp^J \ — -2Y = o / o ^ + a;5,o- (41) 



,2 , p2 
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where the constants x^fl and x^^ can be fixed to by translation invariance in X4 
and X5 (again, the above only represents 1 branch of the solution and must be glued 
on to another branch). This gives simple expressions for L4 and L5 

RpiTT 

La 



^4 



= (42) 
Further, we should note that this gives that 

Ll + Ll = f^. (43) 

In the X4, X5 plane, the ends of the D7 brane have a fixed asymptotic separation 
(not including the warp factor); similar behavior was found in [19] for D5 brane 
embeddings. This gives that L4 and L5 do not parameterize the 2 dimensional class 
of solutions that we have. We can, however, use po = uq = {pi + p|)® and the unit 
vector defining the separation direction L4/L5 = p^/p^ as the two physical quantities. 
Note that one can take fixed L4 and L5 while varying Uq simply by rescaling p^ and 
P5 by the same factor. 

Before turning to a discussion of these embeddings, we note some general points 
about how to find the scaling of Uq with various Lj in the problem. We note that in 
the extremal limit, the action becomes 



Sw = -^-^ / dpu'd {dx^f + {dx^)' + — (44) 

where in the above we have not fixed m = p as a gauge choice. We note that there is 
a Lie point symmetry (for more on Lie point symmetries, see [30]) of the equations 
of motion following from the above action. Namely, if we take 

u C X u (45) 

the action scales as 

5^7 ^ C^Snr- (46) 

This is clearly a symmetry of the solution space of the equations of motion. The 
above symmetry is therefore a local map from solutions to solutions, the definition 
of a Lie point symmetry. 
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One could further scale the embedding coordinates x'^{C,) — C"''x'^{$,) (no sum- 
mation) if one wishes (this rescales the effective V4) and possibly promote the above 
symmetry to a Noether symmetry, however, we will not do this here El- Note that in 
such a scaling, we have not rescaled X4 or 0:5. Therefore, when applying this sym- 
metry to the space of solutions, we know that uq —>■ C x uq however, L4 and L5 do 
not scale at all (as they are related to the asymptotic separation of X4 and X5). This 
further accounts for the same behavior for the D5 embeddings found in [19]. 

One can use such arguments to get the correct behavior in a host of situations. Of 
particular interest is the situation where the metric is of the form gij{u)dx^dx^ +f{u) x 
{du^ + u^dQ\.j) where all metric components along the i,j directions are functions of 
u only, and dQ\.j is the metric of some compact space spanned by coordinates 6^. We 
further assume that the dilaton is a function only of u. The reduced one dimensional 
Lagrangian with = ^* and 6** = 6^{r) will always be of the form G{u) a/(1 + w^dfi^) 
where now all of the ds in the metric dflj^ are to be read as derivatives with respect 
to u. In such a situation, u ^ C x u is a. Lie point symmetry of the action if 
G{C X u) = G{u) X C^, i.e. G{u) = u^. This covers most extremal brane type 
embeddings. Further, such behavior is intuitive for "fiat" type embeddings of the 
probe brane in supersymmetric situations. For a supersymmetric situation, the probe 
brane does not "see" the background it is in, and simply goes along a straight path 
(straight in the sense of the fiat geometry with no backreacted branes). Therefore, the 
angular separation is always a "north pole/south pole" type (at least asymptotically), 
and will not depend on uq at all. As an example, see the discussion of [4] (equation 
(8)). Of course the angular separation may be changed by some constant if there is 
a conical singularity at the origin of the transverse space (at u = 0), as in the more 
recent calculation [31] (see the last paragraph of section 3), however there should 
still be no dependence on uq. 

We can explore such possibilities in general. Consider a situation where the 
metric is of the form Cir"'dx'^dx^'r]ij + G2r"^dr'^ + Gsr"^^'^dQ'^, and the dilaton is of the 
form exp(</)) oc r'^. If one reduces the the DBI action to a one dimensional action via 
X* = x*(r),^* = 9^{r) (with an appropriate number of static gauge coordinates), one 
will always have the Lie point symmetry r — > C x r, — C 2 x*, 6^ ^ 6'^. In such 

m — 71-1-2 

a situation, the asymptotic separation in the is Ax* oc Uq ^ , and the angular 
separation A6^ = constant. For example in the case of the extremal Sakai-Sugimoto 

•^This is because the particular form of this symmetry is dependent on the fact that at least 4 of 
the fields x^,x'^, are along the world volume of the D7 and aligning them to be the first 4 a;^ via 
a Lorentz transformation, and then further gauge fixing these to be . Thus, if we truly wanted to 
extend the above symmetry to a Noether symmetry, we would need to include some combination 
of coordinate transformations and Lorentz symmetry too, which we do not wish to consider here. 
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model [7,8] n = 3/2, m = —3/2 and so Ax^ oc Uq '^R^, where we have used R to 
make the correct units. This answer is indeed correct [16] (see equation (2.18) of this 
work) . 

Our situation is somewhat special, where both the sphere and the flat space 
component have the same warp factor up to a constant. This is the reason that in 
our extremal case, L4 and L5 do not depend at all on uq. It is a statement that a 
pure rescaling in m is a Lie point symmetry, without requiring extra rescahngs of the 
functions x*. Further, in such situations, the boundary conditions are not affected 
by varying uq. Therefore, one may expect such a mode to be normahzable, and 
correspond to a massless mode. One usually thinks of such a mode as a Goldstone 
boson of some spontaneously broken symmetry. Here, the symmetry is the Lie point 
symmetry promoted to a full Noether symmetry, as discussed in the footnote earlier. 
This is a full (non linear) symmetry of the DBl action. It would be interesting to 
see what the symmetry of the DBl action corresponds to on the fleld theory side. 

We now turn back to the embeddings at hand. Note that using the SO (2) sym- 
metry of the background (mixing X4 and X5), one may rotate and translate to a 
frame where X4 = 0, and further, 0:5 (tt) e ' ' ' Such a frame change will 

mix any periodic identiflcations one has made in X4 and X5 (by a slant identiflcation 
"r"), however, the local physics of the embedding will not be dependent on such 
considerations. Therefore, the point X5 = corresponds to where the embedding has 
reached Uq. Therefore, we may summarize the embedding as 




(47) 



This is suggestive of a change of variables. We make the change 



y 



V 




(48) 



(49) 



In these coordinates, the metric becomes 



R 



[rj^j^dx^dx^ + dx^^ + 



9(|/2 + -y2)6 



R 



J {dy^ + dv^) 



+R {y^ + v^^dnl 



(50) 
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We change to "7 brane frame" via 



Gij = exp 9iJ- (51) 

Changing frame in this way gives that det Gp = ^/det gpexp{—(j)) where Gp is the 
puUback of the above G to the D7 world volume. In this frame 

1 1 

= (y'+j')' U^^dx'^dx" + dxl) + -(dy^ + dv^) 

Ri R 9(^2.^^2)8 

+R{y^ + v^)~^ dQl (52) 

The embedding of the D7 is now given by y = ps = Uq, = 0. This gives the 
puUback metric to be 



1 



-ds% = ^ ° „ ^ {r]^,dx''dx'') + {dv^) +R{ul + v^) « dnl (53) 

Ra R 9(^6^^,2)8 

One can now see why y = is a good embedding for any uq: the determinant of the 
above metric does not depend on the warp factors {uq + v'^)\ We will see the solution 
associated with shifting uq in the linearized equations of motion, which we will study 
in the next section. We now add a word of caution. The above symmetry scaling 
Uo to different values is only a classical symmetry. This may not be true quantum 
mechanically, in analogy with the expectations for u\ expressed in section 1 [20]. 



3.1 Scalar Mesons 

We now turn to the fluctuations of the above embedding. To do so, we expand the 
fields as y = Mq + w{v, x^), X4 = + x^^v, x^). Therefore, we find 

^dsQ = — {dx^dxn) H f (dv^) 

i?4 R 9((^3 + y;)2 + ^2)8 ^ 

+R{{ul + wf + v^Y dnl (54) 



R , „ , „^2 . {< + V 



2^| 



- {dywdv + d^wdx^) H {d^x^dv + d^Xidx'^Y 



9(m[S + w2)8 R 
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We wish to expand det{Gp) to second order in w and x^. With the above parametriza- 
tion, this is now easy, as we may think of the third hne as a perturbation of a metric 
given by the first hne. Therefore, we define 

^dsi = i- — idx^dxa) H T (dvn 

+R{{ul + wf + v^Y dnl (55) 

and use this to define a metric Gp. Now, to leading order, the expanded metric 

det(Gp + h) = det(Gp) (l + Gl^'^^K.a,) • (56) 

In the second term of the above, we may substitute u; = into the metric Qp as h is 
aheady second order in w and x^. Further, det((?p) does not depend on w or X4, and 
so there is no need to expand further. We therefore find that, to second order 



/ 



+ r T ,, -dawd^'w + -{d,w)A. (57) 
9(Mg + u2)2 ^ 2^ " ' j ^ ' 

The second order action we wish to consider is therefore 

S = -11^ [ dvd^xi^l-d.wd'^w + lid^wf + Id.x.d'^x, + 9d{d,x,) 
Qs J yyz z z z 



with 



Uq + V 



2 



The above action has positive definite Hamiltonian, so we expect the embedding to 
be stable. 

The eigenvalue problem we wish to satisfy is the following 

dv{9il + V')dvi^,,M.) = -Mli^^^uAV) (59) 
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where the 4D masses are given by mj = ^ and we have defined V = These have 
the general solutions 

/ a/9 - 4M2 1 \ 

+ C,.g(^^-1;-.V) (61) 
,...iV) ^ C,.(l.V^).^.(^.^^,^-^^;l;-^M(^ 



+c2,.\/(i + ) + — — , ^ — , 2^ -y 

where P and Q are Legendre functions, and 2-^1 is the hypergeometric function. 
There are possible bound states in the regime 

< m^,^ = < (63) 

Curiously, we will find this to be true in the next section as well, where we consider 
an embedding into the near extremal background. 

We now change the eigenvalue problem to Schodinger form to analyze the above 
differential equations. We do so by changing coordinates and variables 

V = smh{K) 
{K) = Vcosh(JO^.,Mji^) 

^.mM) = ^-l—- <i^.,MM)- (64) 
A/cosh(ii) 

This brings the equations to the following form 

OI^^^mAK) - Q - l^^^^^ ^^^mAK) = -^Ml^iJ^^MAK) (65) 

d^K^.M.iK) - (1 + ^.MK) = -Im'.^.mM)- (66) 

The second line clearly has no normalizable bound states, where as the first at least 
has the possibility. In fact, the above set of differential equations is exactly the same 
as those encountered in the next section, with a simple modification. 
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In the next section we show that there is only one bound state for the equation 
of the first type above, and in fact this is the the = mode anticipated by the 
discussion above. For = 0, we find 

(r^^ Ci^u, , (^2,,^ sinh(/r) 

*«,,M„=o(i^) = , + , ■ (67) 

ycosh(A) A/cosh(A) 

The first is clearly normalizable, where the second is clearly non normalizable. Fur- 
ther, taking the normalizable solution, we find that the profile in ipw,Mw=o = Ci,io, 
exactly corresponding to shifts in the definition of uq. We will fully analyze the above 
type of differential equations in the next section. 

3.2 Vector Mesons 

We now analyze the fluctuations of the world volume vector gauge field. As in the 
previous section, we find that 

= -m j dQs j dx^dve-'t'^p {l + \Fa,,,Fa,,,gl''''g';'' + ■ ■ ■ ^ . (68) 

which for this case gives 



,3 



1 „ „ 3 



S = -^7^3 / dx^dv- I 1 + l—^l^^-^F^^F^" + 1 --' ^ 0' ^ . p p^.^ ^ 



(69) 

where we now use rj^'^ to raise and lower indices. The equations of motion read 

R'd^'F^, + [ul + 1?) ^ [{ul + v"") ^ = (70) 

= 0. (71) 

As before we expand 

A, = Y,Bl::'\x)Mv), = $^7r„(x)0„(t;) (72) 

n n 

and gauge fix d^^B^ = to find 

J2 R^d^d^B^J'H^ +{ul + v^)-^d,(^ [ul + t 5 ^ = (73) 

n 

({ul + v^) ^ 9,7ro</.o) = ^ (/.o oc ^ ^. (74) 



''0 
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and for 7^ 0, = 0. 

We may change the ipn equations to Schrodinger form via 

V 

= smh.(K) 



^^n){K) = (75) 

(cOsh(K))2 

and taking d^'d^.B^'"'^ = m^E^") = we find 

aim) - (9 + = -ImMK). (76) 

This obviously has no normahzable modes. This means that the only 4D meson 
from this sector is the massless pion ttq. We will see the same phenomena in the next 
section. 



4 A simple antipodal embedding 

In this section we will be concerned with the antipodal embedding P4 = 0, P5 = 0. 
This is the rare exception when we have an analytic solution that comes into the 
region u ~ u\. We will find it convenient to start with the z coordinates introduced 
earlier 

u^ = ul + z\ (77) 

In these variables, the metric reads 



and we make the following coordinate change 



w; = zsin ) (79) 
R 



v = zcoii[^) (80) 



22 



so that 

ds^ = ^""^^ {7]^,dx^dx' + dxfi + R^f^il + ^dnl 

+ ^^={dw^ + dv'). (81) 

Now, the solution to the embedding reads x^i^v) = w{v) = (we will use v as the 
world volume coordinate). 

4.1 Scalar mesons 

We now expand the action to quadratic order in X4,{v,x^),w{v,x'^). This is most 
easily accomplished by going to a "seven brane" frame 

Gu = exp (-^l) gij. (82) 



Changing frame in this way gives that -y/ det Gp = y^det Qp exp(— 0) where Gp is the 
puUback of the above G to the D7 world volume. We first expand the modified line 
element to second order in fluctuating fields x^ and and we find 

gjdsl; = Chi^jydx^dx" + R'^C^dnl + C^dv"^ 

+(^ {dyX4^dv + dfj^x^dx^)^ + C,~i {dywdv + d^wdx^)^ 
3 ^ 5 if ^ , „ , „ ^98 ^ 5 ii?^ ,^0 5 ^ 13 u;^ , o , ^ 



with 



C-^. (84) 

Note that the C, defined in this section is different than the C, defined in the last 
section. 

Next, we expand the action to second order and find 



j d^x^dvdVt^^/^p =^ij d^x^dvC 



1 + ^<9^X4<9^X4 + ]^C,{dvXif 



i5) 
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which has positive definite hamiltonian, and so we expect that this embedding is 
stable. 

The proper eigenvalue problem is 

d. {Cd^^^^mM) = ( - ^) ^^,m^ (86) 

where the 4D masses are rrix = ^ and rriw = The solution to these are the 
following 



i^xM^V) =Cx,i2Fi\--^ '4 + ^^ 4 '2'~ 



2 



/ r p i V^-^M^ i 

^, \/9-4M^ 1 



where 

V = — (89) 

and 2-^1 is the hypergeometric series, and P and Q are Legendre functions. Note 
again that the definition of the above V differs from the last section. 
There are possible bound states in the range 

0<m,^ = ^<^i. (90) 

If or exceed this range of values, the indices of the above functions become 
complex. One can show that around = oo these go as appropriately normalized 
plane wave solutions, and so represent a continuum of states. This corresponds to a 
mass scale at which the scalar mesons may no longer be considered 4 dimensional. 

We expect that the indices of the hypergeometric equation being integers is spe- 
cial, and we find that in the case = that 

'^x,M.=AV) = y2^i ■ (91) 
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In such a case we see that Cx,2 = is a normahzable solution. We will now give an 
argument that this is the only bound state of the above eigenvalue equations. 

We will change the one dimensional problem in hand to Schrodinger form. To do 
so, we take the change of variables and functions 



V = sinh(K) (92) 

^/cosh{K) 
(cosh(K))^ 



^^,M.W = %^^ (93) 
y cosh(A 

ipxMK) = , ,;3 - (94) 



Then, the equations above become 
. /9 1 



dK^ / 4cosh(if)2 

. /9 



*«,,Mji^) - 7 + -:—t:7t^ ^n.,MM<) = -M^^^,mJK) (95) 



,^.^.mM) - - -,^^;^)^.mAK) ^ -M^^.MK). (96) 

In the above problem, we are looking for bound states whose energy is interpreted 
as ii^ = Mf. Clearly there are no bound states for the potential of the first kind (it 
is a hill, rather than a well). In the second case, there is at least the possibility that 
there are bound states. Using the WKB approximation, we take 

cosh" 




2 vrA^^^,-!^ 3 ^ V 4 4cosh(K)2 

V V-12m| + 27 

There is a bound state when n is an integer. Clearly < < by considering 
the maximum and minimum of the potential, and so we plot the above function in 
this range in figure [61 The plot of n crosses an integer value only at = 0. This 
crossing is close to the exact solution of = V^, but not precise. We may expect 
this low level of accuracy for the least energetic bound state, but generically, the 
accuracy of the WKB approximation is improved for higher energy bound states. 
Therefore, we trust that there are no other bound state solutions. 

One further note is in order. For the above configuration, there is no scaling 
symmetry of the DEI action (as there was in the last section). Therefore, it seems 
plausible that the above mode is a pseudo goldstone mode associated with the scaling 
symmetry. Because the symmetry is not exact (due to the function /), the mass is 
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Figure 6: Plot of n(M^). Note that the only integer solution is n = 0, and that this 
approximates the known solution of = \/2. 



lifted from 0. It would be interesting to find a way to track the mode from small to 
large values of Mq/wa and see if this is indeed the case. This would further track from 
being a pure w excitation (when uq 3> ux, as in the last section) to being a pure X4 
excitation (for uq = u\, in this section). Finding the correct linear combination that 
was a bound 4D state, and diagonalized the mass matrix in the general case would 
be instructive. 



4.2 Vector mesons 

As in the previous section, we find that 



(98) 



which for this case gives 

-2 



9s 
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d^F^^ + d.F^'' = ^ a^F^, + ( "^^^^F,, ) = (100) 



The equations of motion read 

a^F^^' = ^ a'^F^„ = 0. (101) 

As before we expand 

^/^ = 5Z4"n^)^n(^^), A = 5^7r„(x)0„,(^) (102) 

n n 

and gauge fix d^B^ = Q to find 

(22 \ 
"^A^d^i^nj Bir^ = (103) 



2 



5. ( ^^^;^5.vro0oJ = ^ 00 oc (104) 



In this case, the eigenvalue problem is 

dv ((1 + V')dvMV)) = -MlMV). (105) 

where again V = —. 

We change to Schrodinger form via 

V = sinh(K) (106) 

MK) = ^^ ^4K) (107) 
y cosh(A) 

and find 

9k-U^ + — ) *n(i^) = -Ml, „vi;„(K) (108) 



4 V cosh(A') 

which obviously has no normalizable states. Note that the eigenvalue problem is 
actually the same as the problem above (with a shift in the potential). Hence, 
the only 4D meson from this sector is again the massless pion ttq. 

4.3 Baryons 

We now turn to the question of determining the size of the baryon, and for concrete- 
ness, we will model this by an instanton in the U (N) gauge field living on the brane. 
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The size of the baryon will be stabilized by the presence of the CS term as in the SS 
model [7,32]. We will only use the above quadratic action plus the CS term arising 
from the term 

Scs = Tf^7 I Tr (exp (F)) A C2 (109) 



where C2 is the potential leading to F3 of the background, and further F = dA — 
i[A^A\ is the U(A'^^)field strength. The appropriate term is the F^ term. We rescale 
F by the canonical factor of 27ra', and integrate by parts to find that 

Scs = TK,j^ {2Tia'f I uj,{A) A F^ = f u;,{A) (110) 



where 



u,{A) = Tr ( AF^ - ^A^^F - ^A^ ) (111) 



and the integration is now five dimensional. This is exactly the term found in [32]. 
We now turn to the quadratic action 



= J dx'^dv^Tr U^F^^ + 2^^^^ ) (112 



and perform the coordinate transformation = Rx'^, v = u^V, and further scale F 
by 27ra' to find 

9s J ^ 

I dxW^Tr {F,^F^^ + 2 (1 + V') F,yF\) 

(113) 



2(27r)3v^ Qs 



where /i, v indices are still raised and lowered by 77^,^. We will drop the " in the 
following, knowing that V is measured in u\ units and x^ is measured in R units. 
Note that in the limit that V ^ Q that the above action goes to that of the free 5D 
Yang-Mills (YM) theory, as it does in [32]. We therefore write down the combined 
action 

S = Sf2 + Scs- (114) 
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Next, we wish to consider the mass of the baryon to leading order in a "small 
size" expansion. We will do this for the SU(2) BPST instanton by checking the 
energy as a function of the instanton size. The BPST instanton configuration is 

F^= = j^^^2'^^^^''^ p^y= ^/fp2y ^- e = V(5^-xo)2 + (v^- WIS) 

The constants {xo, Vq) parameterize the location and p parameterizes the size. Above, 
r* are the Pauli matrices that satisfy Tr(r"r^) = 25"''. We plug this into the La- 
grangian density, and integrate to give the energy of this configuration (measured in 
units of R, as we must "unscale" time to get something with units) 

Eip,V,)R = '"-^ ^ [ d'xdvl-Tr {F,,r^ + 2 (l + V') F,yF\) 

2(27r)'^Vtt 9s J ^ 

Sp^TT^ (2 + V^) 







2(27r)H 




Ua 


{9sN,)^ 


2{2tt)\ 


fa' 9s 


Ua 


i9sN,)l 


2(27r)3\ 


9s 



dV- 



{{v-voy+p'r- 

(47rVo^ + 27rV^ + Svr^) (116) 

which is only minimized at Vq = and p = 0. 

Note that the above energy is the mass of the baryon in this model, which, 
plugging in i? = y/W^xTgl above gives 

This is exactly reproduced by the energy of a D3 wrapping the (sitting at u = u\) 
in the geometry. 



K3 / e-^det(^p) = {2Ti)-'^{a'f I dx^dVls^jR^ 
J J u\gs ^ 



= (27r)-3(a')- V— MA / = ^ / (118) 
9s J 27rQ;' J 

Before we address the effect of the CS term on the size of the baryon, we first 
wish to identify a perturbative paremeter that will make the term parametrically 
larger than the CS term. To do so, we note that the coefficient of the term as 

= ^^^i / '^^'''i^iTr (F,„F- + 2 (1 + V^) F,,F'y) 
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1 



and so the relative coefficient is 2(2^)^ ^ ^- ^^^^ large parameter 

that we will be expanding in. We expect corrections in the size of the instanton to 
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be of order A 2 . With this in mind we scale the spatial dimensions by 

x' = x-hU e {1---3},V = X~W. (119) 



The discussion now exactly mirrors that of the discussion in [32], with the replace- 
ment in their formulas of a = 1/24. Further, the parameter Z in their equations 
must match the appearance of Vq in our calculation above. Therefore, we simply use 
their result with the appropriate changes 

M, = ^(l + \-'(Yl + t + ]+o{\~')] (120) 

where above we have emphasized the scaling with A by including the tilde, i.e. Vq = 
A 2 X Vq, and as we have mentioned before, a = This gives that 

~^ ^ ^ ^ ^ (121) 

Vbn ^/bn A 

To restore units, we note that one must include a factor of R in the above to get the 
order of the radius in the directions. This gives a final result of 

o 6 1 „n A r, R I . . 

pI = ^jR' = (122) 
VStt a V 5 Ma 

Physical length in this background near the point u = u\is measured in dx^ /{R/ua) 
units. Therefore, the physical size of the above distribution is ^ ~ 4.2 in string 
units. This is not very big, and so string corrections may be needed. 



5 Deconfined phase. 

One may also wish to study the physics of the phase where the blackening factor / 
appears on the timelike component of the metric. For such a background, we have 

u R dv? 

ds^ = - (-f{u; UT)dt^ + 6ijdx'dx^ + dxl + dxT) + -— + Rudflj (123) 

R ' u f{u;UT) 

and 

u 1R?' 
exp0 F3 = ^3- (124) 
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In this case ut implies a horizon, however, classically the temperature does not depend 
on this parameter at all (as the periodicity of 0:5 in the last section). Here 

St = 2'kR = ^, (125) 
and further, we can take ut to be a measure of the energy density via 

i^K)=4 (^^)5[l.)4^2 ^ S = 2nRE{ur). (126) 

Above we have used the energy as measured for a near extremal D5 in flat M^'^, and 
S is obtained via the Bekenstein- Hawking formula. 

As in the last section, one can take the embedding of a D7 which is transverse to 
some combination of UjX^x^. Possible configurations are shown graphically in figure 




Figure 7: The configuration of D7 branes in the near extremal D5 brane background 
with the thermal factor on the time component of the metric. The dashed lines show 
a possible configuration with the branes falling into the horizon. 

Again, choosing u to be the world volume coordinate of the brane p, we find the 
reduced action 



-V,Vs [ dpp\ f{p;ur){{dx,y + {dx,y) + ^. (127) 
9s J V P 
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As before, we find the conserved quantities 



^f\p-ur){{dx,f + {dx,Y) + f 

f{p;UT){{dx,Y + {dx,Y) + ^ 



P^ul (128) 
P^u% (129) 



which allows us to solve for 9x4, dx^ as 



dxi = R , ^ (130) 

dx5 = R , (131) 

VpV(p;1)(pV(p;1)-(p^ + p|)) 

where in this section we define 

P = -, d = §, (132) 
Ut op 

This allows us to immediately identify 

where po is the minimal value to which the brane falls (po = Po/ur = uq/ut)- We 
note that the action is a function only of po- The above configurations always give a 
connected brane antibrane for > 1. 

Further, we should note that because the functions only depend on po, 

that for fixed boundary condition at infinity, varying ma does not affect po- This is 
because to maintain the boundary conditions at infinity, one must vary uq in just 
such a way as to keep uq/ut = po fixed. 

One may now ask whether the connected configuration (considered above) is 
preferred or whether the configuration where the brane anti-brane pair simply fall 
into the horizon is preferred. To measure this, we subtract the two (infinite) actions, 
considering the boundary conditions as being placed at some finite (but large) value 
of M, and then relaxing this position to infinity. This defines our regularization 
procedure. 
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The resulting difference we wish to consider is therefore 



dpfP 



PO 



pV(p;1)-p^/(/3o;1) A 



dpp^ 



/(p; i: 



PO '^^1^'VpW;1)-pWo;1) 



-(p^-l).(l34) 



We cannot evaluate this exactly, however, we may plot it numerically. The plot 




Figure 8: The difference in actions between the "U-shaped" and 
urations plotted as a function of po = 



|-shaped" config- 



shows that for ^ < p^ ^ 2.4456 that the separated D7 D7 is preferred, and that 
for > Po ~ 2.4456 that the connected "U-shaped" configuration is preferred. 
Because the temperature is not a function of ut, we may not associate this with a 
phase transition, even in the space of energy density, of which ut is a measure. To 
stress this point, we may easily plug this into the equations for L4 and L5 and find 



(135) 



We note that in the above R does not vary as we vary ut, and neither does L (as this 
is a boundary condition we are imposing). We find therefore that the chiral symmetry 
breaking/restoration is purely a function of the boundary conditions imposed, and 
find that for ^R < L ^ 1.068i? that chiral symmetry is present in the deconfining 
phase, and that for L ^ 1.068i? chiral symmetry is broken in the deconfining phase. 
However, we note that if L < |i? that chiral symmetry is again broken (unless there is 
some problem with an open string tachyon, to which our analysis is insensative) . We 
plot the curve of ^^^^ along with the critical value ^ 1.068, and the asymptotic 
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Figure 9: ^ /?, along with the critical value ~ 1.068, and the asymp- 

totic value of = I ~ 1.047, plotted as red green and yellow respectively. 

Although we have not plotted it, the remainder of the graph for ^^^^ goes to zero as 



value of ~ f ~ 1-047 in figure ([9]) to make this more clear. One can therefore 

see that outside the window < L ^ 1.068i?, chiral symmetry is broken. When 
L is above this window, the brane dips far enough towards the horizon that the 
separated brane configuration is preferred. However, this also occurs below ^R: if 
the asymptotic separation is below this value, there does exist a value of po which 
allows for a connected configuration, however in this case po is always less than /)q, 
and so the separated configuration is preferred. It may be, however, that there are 
no solutions with L < |, and that the disconnected solution is only a local minimum. 
One could argue this by starting with L very large and bringing the branes together. 
It is possible to tunnel from the disconnected to connected solution once L ^ 1.068i?. 
Then as one brings the ends closer together, the value of uq runs off to infinity. In 
fact, from the plot in figure [HI it suggests that the effective potential is running off to 
— oo for smaller and smaller values of L. This suggests that the solutions for L < ^ 
are only classically stable, similar to starting with an "inverted wine bottle" 0^ — 0^, 
and taking = as a solution (when, in fact, it will tunnel out because the potential 
is unbounded from below). We therefore take that there are no solutions for L < ^R; 
for ^R < L ^ 1.068-R, chiral symmetry is broken; and that for L ^ 1.068i? that 
chiral symmetry is broken. This similar to the results of [9] where chiral symmetry 
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was broken for L < 0.97R and restored for L > 0.97R. 

6 Conclusions, and outlook 

The above model offers a new dual description of a model in the same universality 
class as large Nc QCD with flavors, different from that of the SS model [1,7,8]. There 
are certain notable similarities and differences. First, the temperature dependence 
of the theory is quite striking: a thermal "black brane" background is only available 
at a specific temperature T = Further, in the deconfined phase, the boundary 
conditions are functions only of the ratio This has the interesting quality of 
making solutions with the same boundary conditions at infinity easy to find: if one 
scales ut by a factor C, then one must also scale mq by C. This has the effect that 
varying u\ with constant boundary conditions yields no phase transitions: either 
chiral symmetry is always broken or restored, simply given the boundary conditions. 
This allows for the possibility of making such a phase transition simultaneous. Fi- 
nally, we should also note the appearance of a finite height effective potential in the 
radial direction, giving a finite number of 4D mesons, unlike the case of SS. 

There are also some striking similarities. The mechanism for chiral symmetry 
breaking in the confined phase is identical to that of the SS model, and easy to 
visualize. Further, considering the baryons as instantons in the U (Nf) world volume 
gauge field on the D7, an almost identical mechanism for the baryon size being 
stabilized arises, and in both cases the size is of the string scale. 

The biggest drawback of our model is the far to simple spectrum of mesons, re- 
sulting from the finite height and flattening of all effective potentials. One may wish 
to see how one could modify this setup to allow for a larger spectrum of mesons, pos- 
sibly for different values of boundary conditions. More specifically, in the confining 
phase, it would be interesting to see what effect turning on P4 would have. Recall 
that in the anti-podal embedding, the minimum of the effective potential is reached 
when the brane reaches all the way down to u^- Further, recall that introducing a 
non zero P4 does not change this: the brane still reaches the minimal value Uq = u\. 
So, in such a situation, more effective world volume of the brane exists near where 
the bottom of the effective potential exists, possibly widening the bottom of the 
effective potential, and leading to more bound states. Such an analysis, however, 
could be quite involved and need to rely on numerical methods. A further drawback 
is the possible quantum mechanical instability of the system. It would be interesting 
to find solutions to this as well, for example whether the backreaction of the flavor 
branes could help cure this problem. 
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Given the vast number of analysis works of the SS background, there remains 
many possibihties for future work. For example, one may be curious to study the 
meson spectrum in the deconfining phase similar to [10], particularly to see if one finds 
similar patterns to the meson melting found via quasi- normal mode analysis [33] , as 
well as modeling the condensation of various modes as in [14]. Also, including bare 
quark mass will have to follow similar lines as those in the SS model [15,34], and 
studies with such deformations. We have also not addressed any possible transition 
between the deconfined and confined phase. In fact, one should imagine tuning the 
temperature (in the confined phase) all the way up to T = Since this is the 

only temperature that admits two different solutions (with the blackening factor) 
this is the only possible temperature at which a phase transition happens. Note that 
quantum mechanically the region T > does not exist due the Hagedorn instability 
(neglecting the backreaction of the flavor branes). Such a phase transition at T = 
would be driven simply by the presence of the flavors (as the two Euclidean actions 
would be identical for the gravity sector). We would then have to subtract the two 
DBI actions with identical boundary conditions at infinity, but with the blackening 
factor on different components of the metric. Inverting the relation between the 
boundary conditions Li and the Pi is, however, not straightforward. We hope to 
return to these questions in the future. 
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